Observational constraints on generalized Chaplygin gas model 
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The generalized Chaplygin gas model with parameter space a > — 1 is studied in this paper. Some 
reasonable physical constraints are added to justify the use of the larger parameter space. The Type 
la supernova data and age data of some clusters are then used to fit the model. We find that the 
parameters have bimodal distributions. For the generalized Chaplygin gas model, we also find that 
less free parameters fit the data better. The best fit model is the spatially flat model with baryons. 



The best fit parameters are: 
z T = 0.395. 



fimo = 0.044, Wco = — 0.881 and a = 1.57. The transition redshift is 



I. INTRODUCTION 



The mounting evidences suggest that the Universe is 
expanding with acceleration 0,0, 0] ■ This reveals that 
the Universe is dominated by dark energy with negative 
pressure, whose energy density fraction is about 2/3 at 
present. The nature and origin of dark energy is still mys- 
terious, leaving theoretic physicts searching for novel an- 
swers. Althoug a cosmological constant may be responsi- 
ble for the accelerating expansion of the Universe and be 
consistent with current observational data, the extremely 
smallness of the value of the cosmological constant brings 
a big challenge to us. It is also possible that the mecha- 
nism at work is dynamical. A dynnamical candidate to 
provide dark energy may be supplied by a slowly rolling 
scalar field, widely referred as "quintessence field" |4(- 
There are also models of a scalar field with non-canonical 
kinetic term, known as "k-essence" (^J or "tachyonic" 
models . Other alternative dark energy models include 
the holographic dark energy model Q and the extra- 
dimensional motivated models 

Among the cosmological community, the consensus 
about our Universe is the so called "concondance model" : 
95% percent of the Universe is composed of dark compo- 
nents. Recently the Chaplygin gas model with exotic 
equation of state p = —A/p was proposed to explain 
both dark enrgy and dark matter . The Chaplygin gas 
model was later generalzied to the generalized Chaplygin 
gas (GCG) model with equation of state p c = —A/p" 
• The novel feature of this model is that it unifies dark 
energy and dark matter in one model. Sandvik and coau- 
thors claimed that the matter power spectrum essentially 
ruled the GCG model out [Tj- However, their analysis 
does not include the effect of the baryons. Beca and col- 
laborators showed that it is important to include baryons 
in the study of large scale structure and the conclusion 
changed when baryons were included ^2 • ^ should also 
be noted that the results in was based on linear the- 
ory of perurbations and neglected any non-linear effects. 

The GCG model has been extensively studied in the 



literature [13, Hlj . The GCG model was also studied in 
the framework of modified Friedman equation in |15| . In 
this paper, we consider both dust like matter and GCG as 
the source. The dust like matter may be just baryons or 
a portion of dark matter. Unlike other studies, we allow 
the parameter w c q — —A s to be less than —1 and the 
parameter a to be in the region — 1 < a < in addition 
to the usual region a > 0. Some physically reasonable 
conditions are then applied to the model to constrain 
the parameters. Therefore, the model considered here is 
more general in addition to be physical. 



II. GCG MODEL 

In a homogeneous and isotropic universe, the 
Friedmann-Robertson- Walker (FRW) space-time metric 
is 



ds 2 
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For a null geodesic, we have 
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where 



f(n) 



n, 



sinh n , 



-1. 



With both an ordinary pressureless dust matter and 
GCG as sources, the Friedmann equations read 



H 2 



8ttG 



Pc 



a 2 3 
+ 3H(p c + p, 



(fin 



0, 



(3) 
(4) 



where the Hubble parameter H = a/ a, dot means deriva- 
tive with respect to time, p m = p m o(oo/a) 3 is the matter 
energy density, a subscript means the value of the vari- 
able at present time. By using the GCG equation of state 
p c = —A/p", we get the solution to Eq. J3J as 
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Because w c = p c / p c = -A/p* +1 , so A = -w c0 yo" +1 . 
Substitute this expression into Eq. J5J), we get B = (1 + 

w c o)ao (1+Q) PcO +1 - Therefore, Eq. © can be expressed in 
terms of w c o and p c o as 



-w c0 + (1 + w c0 



o \ 3 ( 1 + Q ) 



'(t) 



1 !/(!+«) 



(0) 



It is obvious that the generalized Chaplygin gas behaves 
like the cosmological constant when w c o = — 1 and it be- 
haves like the dust matter when w c q = 0. At early times, 
i.e., the cosmological radius a(t) is small, p c ~ (ao/a) 3 , 
which corresponds to a dust like dominated universe. 
At late times, i.e., the cosmological radius a(t) is large, 
p c ~ constant, which corresponds to a cosmological con- 
stant like dominated universe. Therefore the generalized 
Chaplygin gas interpolates between a dust dominated 
phase in the past and a de-Sitter phase at late times. 
This distinct feature makes the model an intriguing can- 
didate for the unification of dark matter and dark energy. 

The GCG equation of state can be derived from the 
generalized Born-Infeld action ^(| 



a/(l+a) 



(7) 



where = dip/dx^. From the above Lagrangian, we 
can easily get p c = C = —A/p". If we take GCG as 
a quintessence field, i.e., if we take w c o > — 1, then if 
GCG is the only source in a spatially flat universe, the 
potential for GCG is 



V{ct>) 
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+ cosh 2 /( 1+ ")^ 3(1 + a) 
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here we set 8irG = 3. 
By using Eq. |JBJ, we get 



Pc = w c0 p c0 



-w c0 + (1 + w c0 ) 



3(l+a) 



-a/(l+a) 



(9) 

If ao/a -C 1, then to the first order of expansion, p c and 
p c are 



Pc = Pc0(-Wc ) 1/{1+a) 1 



Pc = -pco(-Wco) 1/(1+a) 



1 + Wco /a \ 3 ( 1 + a ) 



w c0 (l + a) \ a 



(10)" 



1 (l + Wco)a /ao\3(i+a) 
w c o(l + a) V a J 



From the above expressions, we see a mixture of a cosmo- 
logical constant with a type of dark energy described by 
a constant equation of state parameter a. So the physical 
meaning of a may be given in this sense. 



Follow Chiba and Nakamura we require the fol- 
lowing physically reasonable conditions: (1) The current 
total energy density is non-negative; (2) The total energy 
density is not increasing; (3) The present sound speed c s 
of the system satisfies < c 2 s < 1 because of causality 
and local stability. The first condition gives fijt < 1, 
where ttk = — fe/(aQ-ffo)- The second condition tells us 
that 1 + <7o — 0^ > 0, where the deceleration parameter 
q = —a/(aH 2 ). The second condition is equivalent to 
the requirement that the effective equation of state for 
the total source w e ff > — 1- The third condition gives 
1 — Q k < jo < 4(1 — fife) + 3<7oi where the jerk parameter 
j = a/(aH 3 ). The deceleration parameter q and the 
jerk parame ter j o are similar to the state finder parame- 
ters used in |18| . 

For the generalized Chaplygin gas model, we get 



jo = 1 - &k 



-aw c o(l 



Wco)QcO, 
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3 3 
-1 + -Q m o + ^fc + 2 (! + w c0 )£l c0 , 



(12) 



(13) 



where fl m {il c } = 8irGp m {p c } / (3H 2 ) and Sl c o = 1 — 
fi TO o — fife. By using the above Eqs. (fT2"|l and JEJ, we 
get the following constraints 

n m0 + (1 + w c o)n c0 > 0, (14) 



aw c o(l + Wco)Q c o < 0, (15) 



fimo + (1 + w c0 )(l + aw c o)£lcO > 0. (16) 

Furthermore, we require that f2 m o > and f2 c o > 0. To 
get accelerated expansion, we also require that qo < 0. 
Therefore, we have one additional constraint 



^mO + (1 + W c o)^cO < g(l - Ofc)- 



(17) 



In the literature, the parameters are usually constrained 
to be w c o > —1 and < a < 1. Some authors also 
considered the possibility of a > 1. From Eqs. I|14I17[) . 
we see that it is possible to get w c o < —1 if a < 0. 
In this paper, we consider the parameter space to be: 
QmO = [0,1], w C Q = [-3,0) an d Wco ^ -1, a = (-1,100]. 
In addition to the constraints (|14I17|) on the parameters, 
we also require that the energy density of GCG given in 
Eq. © is not negative. 

Combining Eqs. J3J), Q and ©, we get the transition 
'redshift zt when the expansion of the Universe under- 
went the transition from deceleration to acceleration 



QcO 



- Wc0 + (l + u )c0 )(l + z T ) 3(1+Q) 



a/(l+a) 



(18) 
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III. SUPERNOVA IA FITTING RESULTS 

In this section, we use the 157 gold sample supernova 
la (SN la) data compiled in Q to fit the model. The 
parameters in the model are determined by minimizing 



E 



[Mobs(^) - K z i)f 



(19) 



where the extinction-corrected distance moduli n(z) — 
5 log 10 (dL(.z)/Mpc) + 25, the redshift z = ao/a — 1, the 
luminosity distance is 



djj = Oo(l + z)r\ = a (l + z)sinn 
l + z . 



QqHq 



Hoy/\£lk 



=smn 



dz' 
dz 1 

W). 



(20) 



the dimensionless Hubble parameter E(z) = H{z)/Hq = 
il m + + fife(l + z) 2 , sinn(a;) is defined as sin(x){x, 
sinh(x)} if fc = 1{0, —1} respectively, and cr.j is the 
total uncertainty in the observation. The nuisance pa- 
rameter Hq is marginalized over with a flat prior as- 
sumption. Since Hq appears linearly in the form of 
5 log 10 Hq in x 2 ' so the marginalization by integrating 
L — exp(— x 2 /2) over all possible values of Hq is equiv- 
alent to finding the value of Hq which minimizes x 2 if 
we also include the suitable integration constant. Be- 
cause we assume a flat prior on Hq, therefore alterna- 
tively we marginalize Hq by minimizing x' 2 = X 2 (y) ~ 
21n(10) j//5-21n[ln(10)^/(27r/XV 1/of )/5] overy, where 
y = 51og 10 i?o- To get the marginalized likelihood of a 
parameter, we marginalize all other parameters by inte- 
grating the probability distribution L = cxp(— x 2 /2) over 
all possible values of the other parameters. 

We first consider the special case w c q = — 1, i.e., the 
LCDM (Lambda cold dark matter) model. The best 
fit to the Riess gold data is Sl m o = 0.31 ± 0.04 with 
X 2 = 176.5. The Akaike information criterion (AIC) is 
X 2 + 2 = 178.5 and the Bayesian information criterion 
(BIC) is x 2 +ln(157) = 181.56 Next we consider the 
spatially flat case Clf. — 0, the best global fit gives that 
ft m0 = 0.073+°; 2 ? 5 , w c0 = -0.97T°;£| and a = 3^ 1a.67 
with x 2 = 173.66. Due to the constraints (|14I17|) . it is 
difficult to find contours for the parameters. The AIC is 
X 2 +2x3 = 179.66 and the BIC is x 2 +31n(157) = 188.83. 
For w c q < —1 and a < 0, We get the local best fit: 
Ono = 0.43tg$, w c q = — 1.40±g;5| and a = — 0.55±g;gl 
with x 2 = 174.12. The AIC is x 2 + 2 x 3 = 180.12 and 
the BIC is x 2 + 31n(157) = 189.29. The probability dis- 
tributions of Q m o, Wco and a are shown in Figs. 11131 It 
is obvious that the distributions have bimodal character- 
istics. Therefore the models p c = —A/p a and p = —Ap a 
fit the supernova data almost equally well. The best fit 
of the full model gives that Q m o — 0.0025, = 0.23, 
w c q = -0.9997 and a = 10.3 with x 2 = 173.46. Appar- 
ently, the full GCG model with curvature term fits worse 



than the spatially flat GCG model and the full model 
tends to be curved LCDM model. 

If we think that the dust like matter source consists 
baryons only, then we can add a prior Q m Q — 0.044 ± 
0.004 20] on the GCG model. In this case, for a spatially 
flat universe, the best fit parameters are £l m o = 0.044, 



WcO 



-0. 



JlfcS and a 



1.57 



+0.1 
-0.94 



with x = 173.95. 



Substituting the best fit parameters into Eq. (|18|) . we 
find the transition redshift is zt — 0.395. The above 
results are summarized in Table [I] Therefore, the current 
SN la data does not favor GCG model over LCDM model. 
Furthermore, more parameters fail to give better fit. 



Model 


x 2 


AIC 


BIC 


LCDM 


176.5 


178.5 


181.56 


GCG1 


173.66 


179.66 


188.83 


GCG2 


174.12 


180.12 


189.29 


GCG3 


173.95 


177.95 


184.06 



TABLE I: The comparison between different models. GCG1 
refers to the globally best fit spatially flat GCG model, GCG2 
refers to the locally best fit spatially fiat GCG model and 
GCG3 refers to the spatially fiat GCG model with the as- 
sumption that the dust like matter source is only baryons. 



Probability distribution of SI 




FIG. 1: The marginalized likelihood distribution of f2 m n for 
the spatially fiat GCG model fitting result to SN la data. 



IV. AGE FITTING RESULTS 

Follow 0| , we define the look back time as 

dx 



■x)E(x) 



(21) 



The current age of the Universe is to = ti(oo) = 14.4±1.4 
Gyr [2l|,|22|. The age of an object i at redshift z is given 
by 



i r ZF 



dx 



x)E{x) 



t L {z F )-t L {z), (22) 
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Probability distribution of w 




over all possible values of df. Alternatively, we marginal- 
ize over df by minimizing \ 2 over df which gives df = 
— Y^i\tL(zi) — tQ hs +t° hs (zi)]/n. Because we already have 
four parameters: f2 m o, fife, w c q and a in the model, we 
use Hq 1 = 9.78k- 1 Gyr with h = 0.72 given by HST 
Key project j23|. The observational data for the age of 
cluster sample is given in [2lJ . We reproduce the data in 

Table nn 



z 


0.1 


0.25 


0.6 


0.7 


0.8 


1.27 


U(z) (Gyr) 


10.65 


8.89 


4.53 


3.93 


3.41 


1.6 


C s -U(z) (Gyr) 


3.75 


5.51 


9.87 


10.47 


10.99 


12.8 



TABLE II: The age data from Ref. |21 



FIG. 2: The marginalized likelihood distribution of w c o for 
the spatially flat GCG model fitting result to SN la data. 



Probability distribution of a 




FIG. 3: The marginalized likelihood distribution of a for the 
spatially flat GCG model fitting result to SN la data. 



where zp is the formation redshift when the object was 
born. From Eq. Ij22|l . we see that 

th{z) = t L (z F ) - ti(z) =t - U(z) - [t - t L (z F )), 

where the delay factor df = to — t^zp) gives the infor- 
mation about the unknown formation redshift zf- As in 
[2l| , we also assume that the delay factor is the same for 
all the objects. The parameters in the GCG model are 
determined by minimizing 



tp — tg 
<7t 



E 



[tL(Zi)-t^ s +tf s ( Zi ) + df} 



of 



(23) 

where t° bs = 14.4 Gyr, a t = 1.4 Gyr and a c = 1 Gyr 
|2l|. The nuisance parameter df is marginalized over 
by integrating the likelihood function L = cxp(— x 2 /2) 



For the spatially flat LCDM model, the best fit result 
is fl m0 = 0.20l[J ^ with x 2 = 1-0. Due to the sparse of 
the data, the result is not as good as that from SN la. 
For the spatially flat GCG model, the best fit results are: 
fimo ~ 0, w c ~ — 1 and a — 6.53. Again, the data does 
not favor GCG model over LCDM model. 



V. CONCLUSIONS 

In this paper, we studied the GCG model which pro- 
vides a phenomenological mechanism of unifying dark 
matter and dark energy. We explored a larger parameter 
space for the GCG model. Instead of studying the usual 
parameter range < a < 1, we extended the parame- 
ter space to be a > —1 with some physical constraints 
on the parameters. We found that the parameters had 
bimodal distributions in general. The constraints from 
Type la SN data and the age data of clusters do not 
favor the GCG model over the simplest LCDM model 
from the standards of the Akaike information criterion 
and the Bayesian information criterion. Therefore, the 
current observations are consistent with both the GCG 
model and the LCDM model. Moreover, the flat model 
fits the observational data better. The only benefit of the 
GCG model is that it may provide a phenomenological 
mechanism of unifying dark matter and dark energy. 
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